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ABSTRACT 

The stress evolution process is taken into account in the linear stability analysis of 
standard thin accretion discs. We find that the growth rate of thermally unstable 
modes can decrease significantly owing to the stress delay, which may help to un- 
derstand the quasi-periodic variability of GRS 1915+105. We also discuss possible 
application of stress evolution to the stability of Shapiro-Lightman-Eardley disc. 



INTRODUCTION 



The standard thin accretion disc (|Shakura fc SunvaevHl973T ) 
has been extensively applied to a variety of astrophysical 
systems, such as active galactic nuclei, X-ray binaries, and 
cataclysmic variable stars. The theory of the standard thin 
disc, which is based on the a stress prescription, predicts 
that the disc will suffer both thermal and viscous instabil- 
ity when the pressure is dominated by radiation pressure. 
However, observations of X-ray binaries show little evidence 
for such instabilities except for GRS 1915+105 (hereafter 
GRS 1915). Since the stability is related to the description 
of stress, the standard a stress prescription was modified 
in some previous works to assure the disc to be stable no 
matter whether gas or radiat i on pre ssure dominates. For in- 
stance, iLightman fc Eardlevl (|l9T4l ) introduced a so-called 
f3 stress prescription, i.e., the stress is proportional to the 
gas pressure rather than the sum of the gas and radiation 
pressure. With the ft stress, the disc will remain stable even 
for the case that radiation pressure dominates over gas pres- 
sure. Furthermore, some other stress prescriptions were in- 
troduced by assuming that the stress is relev ant to both 
the t otal pressure and the gas pressure (e.g., Kato et all 
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sure. Thus, the standard a prescription seems to be more 
acceptable than other prescriptions such as the (3 stress. 



|200S| , p. 108). As pointed out bv lGierliriski. fc Pond l|2004h . 



however, the f3 stress would be in conflict with observations 
as follows. Observations show that the colour temperature 
correction remains stable for a wide range of luminosity, 
whereas discs with /3 stress should have a colour tempera- 
ture correction that changes from ~ 1.8 to ~ 2.7 as the ac- 
cretion rate increases. Conversely, the standard a prescrip- 
tion may predict that the correction should remain stable at 
~ 1.8 for varyi ng accretion rates. Moreover, some recent sim- 
ulations (e.g., iHirose et al.ll2009h showed that the stress is 
roughly proportional to the sum of gas and radiation pres- 
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The slim d isc model was introduced by 
|Xb ramowicz et al.1 (|l988l ). in which it was predicted 
that the thermal instability of the radiation-pressure- 
dominated thin disc will trigger a limit cycle behaviour 
between the thin disc and the slim disc. Such a limit cycle 
was confirmed by time-dependent n umerical calcul ations 
(e.g., Szuszkiewic z fc Miller 2 001; |Li et all l2007h . On 
the other hand, iBelloni et al.l l|l997l ) applied the limit 
cycle theory to interpret the observational quasi-periodic 
variability of GRS 1915. Recently, such limit cycle has also 
been appl i ed to young radio galaxies (e.g., ICzernv et al.l 
120091 ; IWul 120091 ). Comparing the numerical results with 
the observational ones, however, we can find the following 
inconsistency between the theory and the observation. 
The numerical results of Li et al. (2007) showed that 
the duration of the outburst phase thigh is less than 5 
percent of the duration of the quiescent phase ti ow , and the 
luminosity of the outburst phase Lhi g h is around two orders 
of magnitude larger than that of the quiescent phase Li ow ; 
whereas observations of GRS 1915 showed that, thigh is 
comparabl e to ti ow and I/high is only 3-20 tim es larger than 
L low (e.g.. lBelloni et al.lll997l ; IWu et al.ll201~oh . Some efforts 
have been made to improve the theory in order to explain 
obser vations either by som e artifical viscosity prescription 
(e.g.. lNavakshin et al.ll2000l 'l or by additional assumption of 
the energy exchange between the disc and corona (Janiuk 
et al. 2002). 

In the present paper, we will investigate this issue 
through another way by taking the stress evolution process 
into account. Our paper is organized as follows. In Section 
2, we describe the basic assumption for the stress evolution 
process. The time-dependent equations and linear stability 
analysis are presented in Sections 3 and 4, respectively. Con- 
clusions and discussion are made in Section 5. 
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2 MODEL 

The stress in MHD accretion discs can be expressed as (e.g 
Balbus & Papaloizou 1999) 

5B r 8B 



Trip =< p5v r 5Vip 



In 



(1) 



where 8v r (SB r ) and Sv lp (8B v ) are the fluctuating compo- 
nents of the velocity (the magnetic field) in r and <p direc- 
tion, respectively. The angled bracket notation denotes the 
azimuthal average. 

The above formula is a general form and is com- 
plicated for calculation. In the study of accretion discs, 
a more practical formula , i.e., the a stress prescription 
( Sha kura k, Sunvaevlll973l '). is widely adopted: 



-aptot , 



(2) 



where ptot is the total pressure, and a is regarded as a con- 
stant less than unity. With the above simple form, the a 
stress prescription has been extensively applied to the re- 
search of accretion discs, even in the disc stability analysis. 
But in the real system, when the pressure varies, it will take 
a certain period for the stress to adjust. As mentioned in § 1, 
the stability is sensitive to the description of stress, thus the 
stress evolution process may have essential effects on the 
stability. In this paper, we will modify the standard a stress 
prescription by taking the stress delay into account. 

Generally, the growth rate of the stress should be re- 
lated to r°^ p and r TV , where r^ p is the expected stress and 
T rv is the current stress. Thus we have the following general 
function for the evolution of stress: 



dTrip dT rl p , cxp 



Of 



dr 



rip ) i rip J , 



(3) 



where v r is the radial velocity. In a quasi-stationary state, 
the function / should match the following conditions: 
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(4) 



In the present paper, we adopt a simple form as follows: 

^cxp _ „. 

J \'rip : 'rip) — 



(5) 



where t ps is the timescale of stress delay for varying pres- 



sure. We further adopt r®J p 



-aptot according to the a 



prescription, where a is a constant. Combining equations 
(3) and (5), the evolution equation of the stress tensor is 



di 



dt 



+ V r 



-aptot 



dr 



(6) 



3 EQUATIONS 

We first present the time-dependent equations of accretion 
discs for stability analysis. We consider an axisymmetric, 
Keplerian rotating disc (57 = S7k) under cylindrical coordi- 
nates (r, ip, z). The basic vertically integrated equations are 
described as follows (e.g., Kato et al. 2008): 



as la,, , „ 
^ + ^(^0 = 0, 

r or r 1 or 



(7) 



(8) 
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(9) 



(10) 



n = n Ras + n r 



(itoh 3 



(11) 



where E = 2pH is the surface mass density, II = 2ptotH, 
T rip — 2r rv H, and E are the vertically integrated pressure, 
stress, and internal energy, respectively. T is the tempera- 
ture, H is the half-thickness of the disc, fca is Boltzmann 
constant, p is the mean molecular weight, ran is the hy- 
drogen mass, and a is radiation constant. The quantities 
Qvis> Oadv' an( l Qrad represent the viscous heating, advec- 
tive cooling, and radiative cooling rates, respectively. The 
expressions of E, Q^ is , Q^ dv , and Q~^ d are the following: 



E 



3(1-/?) + 
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where f3 is defined as the ratio of the gas pressure to the 
total pressure, i.e., f3 = n gas /n, 7 is the ratio of specific 
heating, and R is the opacity. 

Since the timescale for establishment of vertical equilib- 
rium is the dynamic timescale, which is much shorter than 
the thermal timescale, we modify the stress evolution pro- 
cess (Eq. [6]) as 



OTri, 

dt 



+ V r 



OTrip 

dr 



-all 



(12) 



4 STABILITY ANALYSIS 

Now we have the set of six equations, i.e., Eqs. (7)-(12) for 
linear stability analysis, and the corresponding six physical 
quantities are E, II, v r , H, T, and T Tlp . With subscript 
denoting the quasi-stationary solution, and subscript 1 de- 
noting perturbation quantities, we define the following di- 
mensionless variables: 

El III V r ,l , Hi Tl T r ipl 

o = =— , c = =— , u = — — , h = — — , m — — , / = — — . 

E ^ n rSl' Ho' To' -alio 

With the assumption that all the perturbation quantities 
are proportional to exp(o;i — ikr), the perturbed equations 
corresponding to Eqs. (7)-(12) can be written as follows: 



— a — ikru = , 



— — -u = iakrl— -) —-1 , 
2D. 2 VfT E 



2h = £ - a 



(13) 

(14) 
(15) 



-g(4 - 3/3 + ^ T )a +^h= -~al - Gam + ^aa , (16) 
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t = P{a + m) + (1 - /3)(4m + ft) 



Wins + 1 



t - (Rips) 



dlnllo 



dlnr ojtps + 1 



(17) 
(18) 

- 2 - 2fi(2f2 + 



where k is the epicyclic frequency defined as n 
rd9,/dr). 

Before numerical calculations, we would point out that, 
for long-wavelength unstable modes, the second term on the 
right hand side of Eq. (18) is significantly less than the first 
term, thus the equation can be reduced to 

Ujt pB + 1 

The reason for the above simplification is as follows. Substi- 
tuting equation ()14[) into (|18|l . we will have 



(rfi) 2 E 



2Q 2 iakrHo 



+ (Ot pE 
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91nr a>ip S + 1 



OJtps + 1 



C, (20) 



For standard thin discs, 91nIIo/91nr and k 2 /(2S1 2 ) are 
around unity. Thus, for long-wavelength unstable modes, on 
the left hand side of Eq. (20), the ratio of the second term 
to the first term is 

k 2 (rfi) 2 £ ~ 



(fi* p6 



airJIo 



91nr 



glrJIo 1 
<91nr Lotps + 1 

1 1 



/ 



wt p8 + l k 2 /(2Q 2 ) 
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(21) 



Thus, equation ((20} can be simplified as 



k 2 (r») 2 S 
20? iakrUo 



U)tp S + 1 



which is equivalent to Eq. (19). 

For non-trivial solutions of a, £, u, h, m, and I, we obtain 
the following dispersion relation from equations (|13[l - (|17[ ) 
and (HU, 



where 

Ci = (tttp s 



(22) 



7-6/3 (3 



- 7 Y- 

7-1 2' 7 



-3 7 1-/3 
-1 4-3/3 



tps=0 : thin lines 
=tth: thick lines 




Figure 1. Variation of the growth rate with the radial wave- 
length of perturbation. The thin and thick lines correspond to 
the standard a stress prescription and the stress evolution pro- 



cess with tp 



t t h, respectively. The solid, dashed, and dotted 



lines correspond to /3 = 0.01, 0.1, and 0.4, respectively. 
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When t pB — >• 0, the above dispersion relation is reduced to 
the form in the a stress scenario (e.g., Kato et al. 2008, 
p. 180). 

We numerically solve the dispersion equation Q22p to 
obtain the growth rates of unstable modes. In our calcula- 
tion, we fix a = 0.1 and 7 = 5/3. There exist three modes 



Figure 2. Variation of the growth rate with the timescale of 
stress delay in the long- wavelength limit (kH —¥ 0). 



for each given radial wavelength of perturbations. Among 
the three modes, one is always stable corresponding to a 
negative oj, and is therefore neglected in the following anal- 
ysis. We study the dispersion relation for two cases: (i) with 
stress evolution process (£ ps = i t h); (ii) with the standard 
a stress (t ps = 0). The growth rates of the remaining two 
modes are shown in Figure 1. In this figure, the thick and 
thin lines correspond to the stress evolution process and the 
a stress prescription, respectively. The solid, dashed, and 
dotted lines represent the solutions with /3 = 0.01, 0.1, and 
0.4, respectively, where (3 = 0.4 is the well known critical 
value for the thermal stability with the a prescription. It 
is shown that the growth rate of thermally unstable modes 
(the upper branch) in the stress evolution scenario is sig- 
nificantly smaller than that in the a stress case. On the 
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tth 



(23) 



contrary, the growth rates of the viscously unstable modes 
(the lower branch) of the above two cases are similar. Here, 
the thermal timescale tth is calculated by 

E = 2[3(l-/3)( 7 -l)- / 3] 1 
Q+ ls 3( 7 - 1) an ■ 

Figure 1 indicates that, even though the stress evolution 
process does not change the stability criterion, it may have 
significant influence on the growth rate of the thermally un- 
stable mode. 

For further investigation, we calculate the solutions un- 
der the long-wavelength limit, i.e., kH —¥ 0. In such case, 
the non-zero roots of the dispersion equation (|22|1 are, 



2AL 



-(4^+i)± 



^ + i )2+2 ^ 2 
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where 
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The above equation shows the relationship between the 
growth rate oj and the timescale of stress delay t ps . 

For the case of t ps /tth 3S> 1, Eq. (|24|l is simplified as 



1 2 - 5j0 1 



to 



2t pa 1 + /? 



(25) 



which means that if t ps is sufficiently larger than tth, the 
growing timescale of thermal instability is around t ps . On 
the contrary, for the case of t ps /tth^Cl, Eq. (|24|> is reduced 
to 



2-5/3 



1 

tth 



2At th 1+0 'tth' ^ 

The above growth rate is exactly the same as that in the 
a stress case, in which the growing timescale of thermal 
instability is around tth- 

Figure 2 shows the variation of the growth rate of the 
thermally unstable mode in the long-wavelength limit with 
tp S for P = 0.01, 0.1, and 0.3. The figure clearly illustrates 
that, for the case of t ps > tth, the stress evolution process 
will have essential effects on the thermal instability since the 
growth rate is significantly decreased. 



5 CONCLUSIONS AND DISCUSSION 

In this paper, the stress evolution process, which was ig- 
nored in previous works, is taken into account in the lin- 
ear stability analysis of standard thin accretion discs. We 
show the variation of the growth rate of thermally unstable 
modes with the wavelength of perturbations for two cases: 
with and without stress delay. We find that the growth rate 
with stress delay can be apparently lower than that without 
stress delay. We also make some analytical approximation 
for the growth rate in long- wavelength limit, and present 
the relationship between the specific growth rate and the 
timescale of stress delay. In conclusion, the stress evolution 
process may have essential influence on the thermal instabil- 
ity by significantly decreasing the growth rate, in particular 
for the case in which the timescale of stress delay is compa- 
rable to or even larger than the thermal timescale. 

In a real system the timescale of stress delay remains 



unclear. We would argue that this timescale is possibly com- 
parable to the thermal timescale as follows. In simulations 
with an initially weak toroidal or poloidal magnetic field, the 
magnetic energy may first experience an exponential growth 
during the first few orbits due to the linear instability, and 
then will be followed by the nonlinear evolution. Finally, a 
saturated quasi-steady state phase may form. Totally, it will 
take 15 ~ 20 o rbits for the system to enter a fully turbu- 
lent state (e.g.. lHawlev et al1ll996l ; iFromang &: Papaloizoul 
l2007h . This implies that the evolution timescale of the mag- 
netoturbulence is around 15 ~ 20 times of the dynamical 
timescale td yn , which can be regarded as a duration around 
the thermal timescale tth (e.g., tth ~ tdyn/a and a = 0.1). 
As mentioned in the first section, when the theory of the 
(24) limit cycle behaviour between the standard thin disc and the 
slim disc is applied to the observational quasi-periodic vari- 
ability of GRS 1915, there exists some conflict between the 
theory and the observation on the duration ratio thi g h/tiow 
and the luminosity ratio Lhigh/Llow- In our opinion, the 
stress evolution process may improve the theory to explain 
observations due to the following reasons. When the flow 
suffers thermal instability, the growth rate of M can be sig- 
nificantly decreased by the stress delay. We therefore can 
expect a relatively lower M for the outburst phase, thus a 
corresponding lower Lhigh/Llow- As a consequence, for a cer- 
tain fixed mass accretion rate supply at outer boundary, if 
M at outburst phase drops, the duration of this phase will 
therefore become longer, thus a moderate value of thigh /tiow 
should appear. In other words, the influence of the stress de- 
lay on the growth rate may result in moderate ratios of du- 
ration and luminosity of the outburst phase to the quiescent 
one, which will be more likely to explain the observational 
results than previous calculations. The further investigation 
on this issue may require time-dependent numerical calcu- 
lations including the stress evolution process. 

Moreover, we would point out the possible application 
of the stress evolution process to another geometrically thin 
disc, namely the Shapiro-Lightman-Eardley disc (the SLE 
disc, Shapiro et al. 1976), which was originally introduced as 
the inner region of a thin disc to provide hard X-ray emis- 
sion. However, it is known that the SLE disc is viscously 
stable but suffers thermal instability. If the stress evolution 
process is taken into consideration in the SLE model, and 
the timescale of stress delay is comparable to the thermal 
timescale, we can therefore expect the existence of a quasi- 
stable SLE disc since the decreased growth rate of thermal 
instability may not completely destroy the disc and the vis- 
cous stability may help to suppress the thermal instability 
in the viscous timescale. 
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